Probing Dynamics of Dark Energy with Supernova, Galaxy Clustering and the 
Three- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations 
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Using the Markov chain Monte Carlo (MCMC) method we perform a global analysis constraining 
the dynamics of dark energy in light of the supernova (Riess "Gold" samples), galaxy clustering 
(SDSS 3D power spectra and SDSS lyman-Q forest information) and the latest three-year Wilkinson 
Microwave Anisotropy Probe (WMAP) observations. We have allowed the dark energy equation 
of state to get across —1 and pay particular attention to the effects when incorrectly neglecting 
dark energy perturbations. We find the parameter space of dynamical dark energy is now well 
constrained and neglecting dark energy perturbations will make the parameter space significantly 
smaller. Dynamical dark energy model where the equation of state crosses —1 is mildly favored and 
the standard ACDM model is still a good fit to the current data. 
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Introduction. The recent released three year Wilkin- 
son Microwave Anisotropy Probe observations (WMAP- 
3) P, have made so far the most precise probe of the 
CMB observations. The temperature-temperature corre- 
lation power is now cosmic variance Hmited up to Z ^ 400 
where the glitches on the first peak have now disap- 
peared and the third peak is now detected, which gives 
rise to a better determination on the matter component 
of the universepj. In particular, the direct detection of 
the CMB EE polarization spectrum and better measure- 
ments of TE spectrum have helped a lot in the determi- 
nation on the reionization depth, which is now lower than 
the first year WMAP predictions and with much smaller 
error bars. This has in turn helped to break the degener- 
acy between the slope of theprimordial scalar spectrum 
and the reionization depth'T, Intriguingly now for 
the fittings to the power law ACDM model a Harrison- 
Zel'dovich spectrum is now excluded to ~ Su by WMAP 
alone, which will have profound implications in inflation 
if further confirmed with higher significance level jlj. In 
the fittings to a constant equation of state of dark energy, 
combinations of WMAP with other cosmological con- 
stant are consistent with a cosmological constant except 
for the WMAP -\- SDSS combination, where w < —1 
is favored a bit more than IctU. The measurements of 
the SDSS power spectrum^ is in some sense currently 
the most precise probe of the linear galaxy matter power 
spectrum and will hopefully get significantly improved 
within the coming few years. If the preference of w > —1 
holds with the accumulation on dark energy this will also 
help significantly on our understandings towards dark 
energy. A cosmological constant, which is theoretically 
problematic at present 0, |1] , will not be the source driv- 
ing the current accelerated expansion and a favored can- 
didate would be something like quintessence 0, Q- On 
the other hand, the observations from the Type la Su- 
pernova (SNIa) in some sense make the only direct de- 
tection of dark energy 0, 0, 0, 0, and currently a 
combination of WMAP+SNIa or CMB + SNIa + LSS 



are well consistent with the cosmological constant and 
the preference of a quintessence-like equation state has 
disappeared^. Intriguingly, we are also aware that the 
predictions for the luminosity distance-redshift relation- 
ship from the ACDM model by WMAP only are in no- 
table discrepancies with the " Gold" samples reported by 
Riess et a/jlj. Although the prediction by WMAP is 
consistent with the measurements form the Supernova 
Legacy Survey (SNLS)^31j we are aware that the 71 high 
redshift type la supernova alone are too weak for a cos- 
mological probe and even when combined with the 44 
nearby SNIa their constraints on dark energy are not yet 
comparable with the Riess " Gold" samplejl^ . Al- 
though the discrepancy might be some systematical un- 
certainties in the Riess "Gold" sample, this needs to be 
confronted with the accumulation of the 5-year SNLS ob- 
servations and the ongoing SNIa projects like the Super- 
nova Cosmology Project (SCP) and from the Supernova 
Search Team (SST). Alternatively, this might be due to 
the implications of dynamical dark energy. 

The rolling scalar field of quintessence typically pre- 
dicts some evolutions of dark energy equation of state 
and with the accumulation of the various kinds of obser- 
vations we are now also able to constrain the dynam- 
ics of dark energy over a constant equation of state. 
Also the WMAP team has considered the constraints on 
the parameter space where the equation of state is less 
than —1. Such an equation of state is originally moti- 
vated by the scalar field model of phantom^^, where 
there is a negative kinetic term and typically w is not a 
constant either for the case of phantom and it is natu- 
ral that we can expect a nonzero dw/dz with phantom. 
The model of phantom has been proposed in history due 
to the mild preference for a constant w < —1 by the 
observations fisj . Although phantom violates the weak 
energy condition(WEC) and faces the dilemma of quan- 
tum instabilities 16], it remains pos sible in the descrip- 
tion of the nature of dark energyjlTj. 

Previously with the accumulation of SNIa data 
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and especially after the SNIa observations from the 
HST/GOODS programjij, many groups have started 
to probe the time dependence of dark energy equation 
of state [H El I23, Hi llil. Intriguingly an equation 
of state which crosses the cosmological constant bound- 
ary are somewhat favored 0, 0, I2(tI l2ll |. Such a kind 
of behavior is nontrivial physically since conventional 
quintessence and phantom cannot realize such a behav- 
ior. In Ref.[23| we dubbed the new kind of dark energy 
quintom in the sense that it resembles the combined be- 
havior of quintessence and phantom. Although the model 
of k-essence which has non-canonical kinetic term^,'2^ 
can both have a quintessence and phantom-like behav- 
ior, as shown in Refs.[2^ |23| a crossing behavior is 
not viable. Mathematically, the crossing reads that there 
exists at least one pivot redshift, namely z* satisfies: 

^\izl)^0,w{zt)^-l (1) 

where z denotes the redshift and w{z) is the functional 
form of EOS evolution. Interestingly, the quintom mod- 
els differ from the quintessence or phantom in evolution 
and the determination of the fate of universe 28j. There 
exist lots of interests in the literature presently in build- 
ing of quintom-like models. For example, with minimally 
coupled to gravity a simple realization of quintom sce- 
nario is a model with the double fields of quintessence 
and phantom^, 2S,i 2^ 1^ <^-> JMJ ■ In such 

cases quintom would typically encounter the problem of 
quantum instability inherited from the phantom com- 
ponent. However in the case of the single scalar field 
model of quintom, Ref. 36] considered a Lagrangian with 
a high derivative term adding to the kinetic energy and 
its energy-momentum tensor is equivalent to the two-field 
quintom model. Such a model is theoretically possible to 
resolve the problem of quantum instabilities and needs 
further investigations [s^ . 

Given the lack of theoretical understandings, the cos- 
mological observations play a crucial role to study dark 
energy. Probing the dynamics of dark energy is of great 
significance to shed light on theory. SNIa observations 
which measure the luminosity distance depending on 
Hubble parameter H{z), are relatively sensitive to the 
dynamics of dark energy. Moreover, dark energy also 
leaves imprints on Cosmic Microwave Background(CMB) 
through the distance to the last scattering surface es- 
pecially when we take the perturbation of dark energy 
into consideration. The fluctuations of dark energy can 
lessen the ISW effect and lowers the power spectrum at 
small multipoles. Besides affecting H{a), dark energy 
also modifies the growth rate of structures via pertur- 
bation equations. Thus the LSS data such as SDSS 3D 
power spectrum and Lyman-a forest, can further be used 
in the determination of the dark energy parameter as well 
as breaking the degeneracy among the various cosmolog- 
ical parameters. 

Given the fact that due to the problems on dark en- 
ergy perturbations '26!| in the conventional parametriza- 



tions of dark energy equation of state one cannot make 
global fittings with w getting across the cosmological con- 
stant boundary. Previously for example Ref. jSg studied 
the perturbations of dynamical dark energy only for the 
regime where w > —1 and Ref. Is^l considered both the 
cases for w > — 1 and w < —1, but did not include the 
perturbations for quintom like dark energy; while other 
global analysis like Ref. |40l l4ll did not consider dark en- 
ergy perturbations. In Ref. |26j| we gave a method which 
for the first time allows to study quintom-like dark energy 
perturbations with parametrization of w, which resem- 
bles the behavior of the simplest double-field quintom^. 
Our method has allowed a global analysis on the dynam- 
ical dark energy equation of state with the observations 
and like the case for a constant equation of statef43, 143, 
which has also been recently shown by Spergel et ij, 
we find incorrectly neglecting dark energy perturbations 
typically make the parameter space smaller and hence 
with notable bias. 

The aim of current paper is to study the current up-to- 
date observational constraints on dynamical dark energy. 
We extend our previous work of Ref. 26] and study the 
full observational constraints on dynamical dark energy. 
In particular we pay great attention to the effects of dark 
energy perturbations when EOS crosses -1 as in (P). Per- 
turbations of the quintom-like models have been studied 
extensively in Ref. (26] . Our paper is structured as fol- 
lows: in Section II we describe the method and the data; 
in Section HI we present our results on the determination 
of cosmological parameters with (WMAP-3) 2^ , SNIa 
[11 m , Sloan Digital Sky Survey 3D power spectrum 
(SDSS-gal) and SDSS Lyman-a forrest data (SDSS- 
lya) 13] by global fittings using the Markov chain Monte 
Carlo (MCMC) techniques i^ESIH; conclusions and 
discussions in are presented in the last section. 

Method and data. In this section we firstly present the 
general formulae of the dark energy perturbations in the 
full parameter space of w{z) especially for the crossing 
modeltH] In the global fitting to Cm'b,SN la and LSS 
data, we adopt the parametrization of EOS as follows: 485 

w{z) = wo + wi— — (2) 
1 + z 

Despite our ignorance of the nature of dark energy, it 
is more natural to consider the DE fluctuation whether 
DE is regarded as scalar field or fluid rather than simply 
switching it off. The conservation law of energy reads: 



Although the multi-field dark energy models are more challeng- 
ing on theoretical aspects of naturalness, given that we know very 
little on the nature of dark energy, the energy momentum of such 
models can be identified with single field scalar dark energy with 
high derivative kinetic terms IStl . Our phenomenological formula 
of perturbations on DE corresponds to such models of multi-field 
(qui ntom ) with a negligible difference around the crossing point 
of -im. 
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FIG. 1: The lines wo = —1 and wo + wi = —1 divide the 
itJQ — parameter space into four parts. The gray-shaded re- 
gions are for crossing models where the EOS of "Quintom A" 
models is greater than -1 in the past and smaller than -1 to- 
day while the "Quintom B" models cross -1 in the opposite 
direction. 



where T'"' is the energy-momentum tensor of dark energy 
and " ;" denotes the covariant differentiation. Working in 
the conformal Newtonian gauge, Equation© leads to the 
perturbation equations of dark energy as follows : 

h = -(l + w)(6i-3$)-3H((5p/(5/9-u;)(5 , (4) 

B = -n{i-'dw)e-^^e + e{^-^5 + ^) ,(5) 

w 1 -\- w 

Sp = c^Sp + TSS . (6) 

where is the adiabatic sound speed and TSS is the 
entropy perturbations. For simplicity we neglect the en- 
tropy perturbations and assume the sound speed to be 
Cg — 126\. For the models where the EOS doesn't cross 
-1, the lower left and the upper right region of FIG^ 
the above equationl^)), ||SJ is well defined. For the cross- 
ing models in Eq.lQJ, graphically the gray-shaded area 
of FIG^ the perturbation equation I^J, ©is seem- 
ingly divergent. However basing on the realistic two- 
field-quintom model as well as the single field case with 
a high derivative term |2^, the perturbation of DE is 
shown to be continuous when the EOS gets across -1, 
thus we introduce a small positive parameter ^ to divide 
the full range of the allowed value of the EOS w into 
three parts: 1) w > -1 + ^; 2) -1 + ^ > w > -1 - ^; and 
3) w <-!-£,. 

For the regions 1) and 3) the perturbation is well de- 
fined by solving Eqs.Q), |(SJ| as shown above. For the 
case 2), the perturbation of energy density S and diver- 
gence of velocity, 9, and the time derivatives of S and 9 
are finite and continuous for the realistic quintom dark 
energy models. However for the perturbations with the 



above parametrizations clearly there exists some diver- 
gence. To eliminate the divergence typically one needs 
to base on the multi-component DE models which re- 
sult in the non-practical parameter-doubling. A simple 
way out is to match the perturbation in r egion 2) to the 
regions 1) and 3) at the boundary and set fM l2fil l5?il | 

S^O , 9 = 0. (7) 

We have numerically checked the error in the range 
\Aw = ^1 < 10^5 and found it less than 0.001% to the 
exact multi-field quintom model. Therefore our matching 
strategy is a perfect approximation to calculate the per- 
turbation consistently for crossing models 1^. For more 
details of this method we refer the readers to our previous 
companion papers 0, |2^ . 

We have modified the publicly available Markov Chain 
Monte Carlo package camb/cosmomc|5lj ^ to allow for the 
inclusion of dark energy perturbations with EOS get- 
ting across -l|2^ and then sampled from the following 
8 dimensional cosmological parameter space using the 
Metropolis algorithm : 

p = {ujb, ^^c,Qs,T, Wo, wi, Us, \og[lO'^° As]) (8) 

where ujb = flbh^ and tOc = flch^ are the physical baryon 
and cold dark matter densities relative to critical den- 
sity, 9s is the ratio (multiphed by 100) of the sound 
horizon and angular diameter distance, r is the optical 
depth, As is defined as the amplitude of the primor- 
dial scalar power spectrum and Ug measures the spec- 
tral index. We have also marginalized over the bias 
factors b defined as 6 = [Pgaiaxy{k)/ P{k)]^^'^ which are 
assumed to be constant. Basing on the Bayesian anal- 
ysis, we vary the above parameters fitting to the ob- 
servational data with the MCMC method. Through- 
out we assume a flat universe and take the weak pri- 
ors as: T < 0.8,0.5 < Ug < 1.5,-3 < wq < 3, — 5 < 
wi < 5, 0.5 < 65 < 10, a cosmic age tophat prior as 
10 Gyr< to <20 Gyr. In addition, we make use of the 
Hubble Space Telescope (HST) measurement of the Hub- 
ble parameter Ho = lOO/i km s^^Mpc^^ 52j by multi- 
plying the likelihood by a Gaussian likelihood function 
peaked at around h = 0.72 with a standard deviation 
a = 0.08. We impose a Gaussian prior on the baryon 
and density f^b/i^ ^ g 922 ± 0.002 (1 a) from Big Bang 
nucleosynthesis ^3 . 

In our calculations we have taken the total likelihood 
to be the products of the separate likelihoods of CMB, 
SNIa and LSS. Alternatively defining = — 21og>C, we 
get 

Xtotal — XCMB + XsNIa + XlSS ■ 

In the computation of CMB we have included the three- 
year temperature and polarization data with the routine 



^ http:/ /cosmologist.info/cosmomc 
yhttp: / / camb.infO| 
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for computing the likelihood supplied by the WMAP 
team P, In the calculation of the likelihood from 
SNIa, we use the 157 "gold" set of SNIa pubHshed by 
Riess et al in Ref.|llj| and marginalize over the nuisance 
parameter. For LSS, we firstly use the code of CAMB to 
generate the theoretical 3D matter power spectra of ev- 
ery model and fit to SDSS 3D power spectra data*^ using 
the likelihood code developed in Ref. [s^l- To compute 
the lyman-Q; forrest likelihood, we use the SDSS lyman-a 
data and corresponding likelihood code^^l ■ 

Results. In this section we present our results and fo- 
cus mainly on the dark energy parameters. In particu- 
lar to show the effects of dark energy perturbations we 
present the resulting constraints on the parameters for 
two cases simultaneously: one with and the other (incor- 
rectly) without dark energy perturbations. 

In Table 1 we list the mean and 1, 2cr constraints on 
dark energy, inflation and rcionization related parameters 
with/without DE perturbations. In our notations in the 
error bars of the 2a constraints we have included the la 
contributions of uncertainties. By virtue of the combined 
data and especially for the new precision WMAP-3 data, 
we find that all these parameters are well determined in 
our both cases. The center value of dark energy param- 
eters illustrate the favored dynamical DE models with 
EOS crossing -1. In our framework of dynamical DE, a 
lower optical depth t is still favored which is consistent 
with the new results of WMAP teamUQ. Moreover, a 
simple scale-invariant primordial spectrum is disfavored 
at slightly larger than 2-a even in the presence of dynam- 
ical dark energy where the perturbations are included. In 
fact from Fig[5]we find r is still well constrained in the 
present case, which breaks the ns — t degeneracy and 
favors nontrivially a red primordial scalar spectrum. 
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FIG. 2: 1-D constrains on individual parameters using 
WMAP-3+157 "gold" SNIa+SDSS-gal+SDSS-lya with the 
model of dynamical dark energy discussed in the text. Red 
SoIid(Blue dashed) curves illustrate the marginalized distri- 
bution of each parameter with/without DE perturbation. 
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FIG. 3: 68% and 95% constraints of the 2-D contours among 
dark energy and the background parameters. Red solid and 
blue dashed lines are for the perturbed and (incorrectly) un- 
perturbed DE respectively. 



In Fig|21 we delineate the corresponding posterior one 
dimensional marginalized distributions of wq, wi, Us and 
T from our MCMC results. Since the parameters, espe- 
cially for wi, are not perfectly Gaussian distributed, the 
mean and the marginalized results for each parameter 
may be different. For the compliment of table (1), we 
find the best fit values constrained by the full dataset is 
Wo -1.363, wi = 1.325, = 0.964, t = 0.082 with DE 
perturbation and wq = —1.212, wi — 0.839, = 0.953, 
T = 0.064 with when dark energy perturbations are (in- 
correctly) switched off. One can find that almost all of 
the best fit parameters have been modified by the effect 
of DE perturbations. Moreover, the allowed parameter 
space has been changed a lot and the constraints on the 
background parameters have been less stringent when in- 
cluding the dark energy perturbations. This can also be 
clearly seen from the two dimensional contour plots in 
FigEl The reason is not difficult to explain. The ISW ef- 
fects of the dynamical dark energy boosts the large scale 
power spectrum of CMB_26]. For a constant equation of 
state Ref. has shown that when the perturbations 
of dark energy have been neglected incorrectly, a sup- 
pressed ISW will be resulted for quintessence-like dark 
energy and on the contrary, an enhanced ISW is led to 
by phantom-like dark energy. In this sense if we neglect 
dark energy contributions, there will be less degeneracy 
in the determination of dark energy as well as the rel- 
ative cosmological parameters. However, dark energy 
perturbations are anti-correlated with the source of mat- 
ter perturbations and this will lead to a compensation 
on the ISW effects, which result in a large parameter 
degeneracy^2- In fact as we have shown that crossing 
over the cosmological constant boundary would not lead 
to distinctive effects'^, hence the effects of our smooth 
parametrization of EOS on CMB can also be somewhat 
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TABLE 1. Mean and I, 2a constrains on dark energy, spectral index and optical depth parameters using the combined data of WMAP-3, 
"Gold" SNIa sample, SDSS 3D power spectra and SDSS lyman-a forest information with/without DE perturbation. Note that in the 
error bars of the 2<j constraints we have included the la contributions of uncertainties. 
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identified "with a constant effective equation of state [sj 
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liowever tfie SNIa and LSS observations wiil break sucli 
a degeneracy -witfi additional geometrical constraints. 
Thus for the realistic cases of including dark energy per- 
turbations, the correlations between the dark energy and 
the background parameters as well as the auto correla- 
tions of the background cosmological parameters have 
been enlarged, as can be seen from FigEl 




FIG. 4: Contour plot of wq and wi. Here we use the same 
parametrization and data sets as FIG. 2, red solid and blue 
dashed lines are for perturbed and unperturbed DE respec- 
tively. 

Dark energy perturbation affect dark energy param- 
eters most directly and significantly, which can also 
be seen from Fig0| on the constrains in the {wq , wi ) 
plane. For the parameters {wq, wi) the inclusion of the 
dark energy perturbation change its best fit values from 
(-1.212, 0.839) to (-1.363, 1.325). Dark energy per- 
turbation introduces more degeneracy between wq and 
Wl thus enlarges the contours a lot. From the figure we 
can see that dynamical dark energy with the four types 
are all allowed by the current observations and Quintom 
A seems to cover the largest area in the 2-dimensional 
contours with all the data we used. 

As shown in FigOl that wq and wi are in strong cor- 
relations. The constraints on w{z) are perhaps rel- 
atively model independent, as suggested by Ref . [4ll| . 



Following we obtain the constraints on w{z) by com- 
puting the median and 1, 2a intervals at all redshifts up 
to z = 2. In FigElwe plot the behavior of the dark energy 
EOS as a function of redshift z, we find that at redshift 
z ^ 0.3 the constraint on the EOS is relatively very strin- 
gent. One can see that the perturbation reinforces the 
trend of DE to cross -1 at z ~ 0.3. However despite in 
some sense w{z) is well constrained by the current data, 
the quintom scenario is only favored at ~ la by the full 
dataset from CMB, LSS and SNIa. The accumulation 
of the observational data are still urged especially when 
for the probe of the dynamical dark energy. We find the 
value at z = 0.3 is restricted at 
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0.3) = -1.003^-.-ro1^^ (11) 
for the case without dark energy perturbations and 

w{z = 0.3) = -imst^iT-'oilt (12) 

when including dark energy perturbations. Correspond- 
ingly at redshift z ^ 1 the constraints turn out to be 



1) = -0.869 



+0.149+0.230 
-0.160-0.407 



without perturbations and 

w{z = 1) = -0.846 



+0.209+0.277 
-0.221-0.789 



(13) 



(14) 



when including dark energy perturbations. One should 
bear in mind that such a constraint is not really model 
independent, as shown in Refs. [58ll59l|. 

Discussion and conclusion. In this paper we have per- 
formed a first analysis on dynamical dark energy from the 
latest WMAP three year as well as the SN la and LSS 
information. Our results show that when we include the 
perturbations of dark energy, the current observations al- 
low for a large variation in the EOS of dark energy with 
respect to redshift. A dynamical dark energy with the 
EOS getting across —1 is favored at around la with the 
combined constraints from the latest WMAP , SDSS and 
the "gold" dataset of SNIa by Riess et al. 

Compared with our pre- WMAP-3 results^ we find 
now the constraints on dark energy parameter space are 
improved significantly. This lies on the fact that the 
CMB observations have been improved significantly and 
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FIG. 5: Constrains on w(z) using WMAP 4- 157 "gold" SNIa 
data + SDSS with/without DE perturbation. Median (central 
line), 68%(inner, dark grey) and 95%(outer, light grey) inter- 
vals of w(z) using 2 parameter expansion of the EOS in (4). 



also an inclusion of the smaller scale power of the Lyman 
a also help a lot to break the degeneracies. We also note 
in addition to the discrepancy between WMAP-3 and 
the Riess sample of "Gold" SNIa, WMAP-3 also prefers 
a lower value of erg. erg is an important quantity charac- 
terizing the size of fluctuations on the galactic scales and 
this will have some profound implications on the studies 
of structure formations. On the other hand Lyman a 
data prefers a higher value of erg and in some sense, our 
current method of global fittings are not very strong in 
the probe of such discrepancies unless all the notewor- 
thy inconsistencies are all due to neglecting the dynam- 
ics of dark energy. While the current WMAP-3 data and 
the WMAP + SDSS combination favor both a deviation 
from ns = 1 and w — —1, where some nontrivial dynam- 
ics might be available "simultaneously" the both phases 
of accelerated expansions more parameters characteriz- 



ing the dynamics and inherent physical quantities like 
the tensor to scalar ratio r of the primordial spectrum, 
running of the scalar spectral index dns /d In k and also 
the running of w as shown in the current paper need also 
to be considered for both probes of the dynamics and also 
to avoid some possible bias due to simply assuming no dy- 
namics in one sector 'SOj. Moreover in our current study 
we have neglected the secondary Sunyaev-Zeldovich(SZ) 
effects in CMB calculations, which can also give rise to 
some minor shifts on the parameter space. But as indi- 
cated in Ref. 0| such shifts are typically small and the er- 
ror bars on cosmological parameter estimations are sim- 
ilar when compared with cases neglecting the SZ effects. 
Further detailed analysis on the implications of dynam- 
ical dark energy in light of the current observations are 
currently in progress. 

In our all results, the perturbation of dark energy plays 
a significant role in the determination of cosmological pa- 
rameters. Neglecting the contributions of dark energy 
perturbation will lead to biased results as shown explic- 
itly. In the next decade, many ongoing projects in the 
precise determination of cosmological parameters will be 
available. We can hopefully detect the signatures of dy- 
namical dark energy like quintom through global fittings 
to the observations, where it is crucial for us to include 
the contributions of dark energy perturbations. 
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